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t— I . Abstract 

' We present a possible scaling solution to pre-equilibrium evolution 

which interpolates between the instability present in the dense gluon 
\^ | system produced immediately after a heavy ion collision and the final 

t— I ' equilibration which occurs later. Our solution depends on a single 

parameter 5. Depending on the value of 5, our proposed solution 
l/^ , matches onto the bottom-up picture either at an intermediate stage 

| or toward the end of the evolution given by bottom-up. We discuss in 

■ detail the reasons why we believe our solution is self-consistent, and 

we also point out why it is difficult to actually prove consistency. 

a.: 
^ ■ 
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1 Introduction 

The bottom-up picture [1] of thermalization in heavy ion reactions was an 
attempt to describe the various stages that the dense system of gluons pro- 
duced in a heavy ion reaction go through on the way to equilibration. The 
picture, which relied on the observation that inelastic processes were im- 
portant for kinetic equilibration [2J, was relatively simple with three distinct 
stages. An early stage dominated by hard gluons having momenta on the or- 
der of Q s , the saturation momentum of the colliding nuclei; an intermediate 
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stage where hard gluons dominate the total number of gluons but softer glu- 
ons dominate the Debye mass calculation; a final stage where the soft gluons 
are thermalized and receive energy, by QCD branching, from the hard gluons 
until all the energy is removed from the hard gluons and the whole system 
is equilibrated. 

Arnold, Lenaghan and Moore[3] pointed out a severe problem with the 
bottom-up picture. Rather than having a screening mass at early times the 
system of gluons produced has a very asymmetric momentum distribution 
which leads to a mass characterizing an instability [3-7]. Thus, the early 
stage of the bottom-up picture is not tenable. The instability clearly leads 
to an increase in the transfer of energy from hard to soft gluons, and so 
one could hope that it produces a rapid thermalization. However, Arnold 
and Lenaghan [8] showed that complete equilibration cannot occur before 
Q s r = oT 1 l z a time not much shorter than that which occurs in the bottom- 
up picture. Thus the instability cannot lead directly to equilibration since 
that would give an equilibration time parametrically on the order of Q s t ~ 1. 
Thus another mechanism might be necessary to fill this gap. There are no lack 
of proposals in the literature for this ranging from a strongly coupled quark- 
gluon plasma[9,10], a mechanism using the Hawking-Unruh effect [11] and 
collinear enhancement within the bottom-up picture[12]. No matter what 
the mechanism really is as long as it is perturbative in nature 4 , it would be 
natural to look for a scaling solution which begins when the instability runs 
its course and ends when equilibration occurs. And, of course, this scaling 
solution cannot be the bottom-up solution. 

In this paper we search for solutions to the Boltzmann equation where 
the collision term is given in terms of two to two and two to three gluonic 
processes. We allow both the hard gluons, the remnants of the gluons pro- 
duced in the initial ion-ion collision, and the soft gluons, produced later on 
in the evolution of the system, to be part of the collision term. As in Ref.l 
we treat the collision term in a dimensional way without trying to keep track 
of constant and logarithmic prefactors to the dominant behaviors. 

We have found a one-parameter family of scaling solutions parametrized 
by a positive number 8. We believe that these are the only scaling solutions 
which exist. They are exhibited in Eq.(l) for N s the number density of soft 
gluons, k s the soft gluon momentum, f s the soft gluon occupation number 
and mo the Debye screening mass. Unfortunately, we cannot prove that mn 

4 Clearly not all the above mechanisms are perturbative. 
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is actually a screening mass rather than an instability mass. We show in Sec. 
2 that our proposed solution only has internal consistency if mo is interpreted 
as a screening mass. As we believe the system should have a scaling solution 
connecting the instability to equilibration, and we believe the solutions in 
(1) are the only possible scaling solutions, we are hopeful, but not absolutely 
sure, that the mo in (1) is a genuine screening mass. At the end of Sec. 2 we 
discuss the competing mechanisms which occur in trying to decide whether 
run is screening or whether it corresponds to an instability. 

In Sec. 3 we show that our proposed solution matches onto the interme- 
diate stage of the bottom- up picture for < 5 < 1/3, with 5 = being the 
bottom-up scenario and 5 = 1/3 having the matching occur just as the soft 
particles are beginning to equilibrate. Solutions with S > 1/3 create equi- 
libration of soft particles earlier than the corresponding time in bottom-up, 
and in this case the matching with bottom-up only happens as the whole 
system is reaching equilibration. Although the thermalization we find here 
is parametrically the same as in bottom-up the system may well become 
"effectively thermalized" at an earlier time as suggested in Ref. [6]. 

As in bottom- up equilibration occurs when Q s r ~ a~ 13 ^ 5 , independent of 
5. The fact that we have a family of solutions rather than a single solution 
reflects the weakness of our approach. An approach utilizing more dynamics, 
such as an actual solution to the Boltzmann equation, should find a unique 
solution matching onto the instability at early times. 

Finally, in Sec. 4 we give a simple proposal of when and how the insta- 
bility may run its course. 



2 The scaling solution 

Our procedure here is just to present a solution for the evolution of the 
pre-equilibrium gluon system, after the instability has run its course, and 
then to check that this solution appears to satisfy all the requirements that 
such a solution should satisfy. The part of the following discussion which 
is somewhat unsatisfactory is our inability to prove that what we call the 
Debye mass, mo, actually obeys m 2 D > 0. At the end of our discussion in 
this section we comment as to why this is not so easy to show. 
Our proposed solution is given by 

AT ^ t & 



«(Q S r) 4 /3-<5' s (Q sT ) 1/3-25/5' 
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afs ~ (Qst) V 3 + 5 / 5 ' mD ~ (Q s r) V2-M/10 W 

with 5 a positive real number and where, as in Ref. [1], the symbol ~ means 
that we are unable to evaluate the constant factors (perhaps logarithmic 
factors) which accompany the dimensional factors, the powers of a and the 
powers of Q s t which we have evaluated. Again, as in Ref. [1], we separate 
hard gluons having momentum on the order of Q s from the "soft" particles 
produced around time r which have momentum k s . In addition, there are of 
course gluons which have been produced in the interval 1/Q S <C t <C r and 
these gluons, produced at r , have 

N s (r,r ) ~ m/w ^( r °) 



a(Q s r)(Q s T )V3-5' (Q a7 . )i/3-M/5> 

f , , (g^o) 1/3+3/5 (0 , 
a M r - r o)~(^jW ( 2 ) 

where, for example, N s (r, r ) stands for particles produced around time r 
and measured at time r > r . In evaluating af s (r, r ) we have used 

faij , To) ~ ;2 /V/ f\ 

because the soft particles produced at r , just as the hard particles, have a 
longitudinal momentum determined by multiple scattering at time r. Since 
N s (r,To)f s (r,T ) < N s f s the particles produced at r are not so important 
for scattering which occurs at time r, and so for the present we neglect them. 
At the end of this section we will come back to the role they play, and we 
shall see that the separation of N s from N s (t, t ) is subtle. 
The solution (1) has the following properties: 



< ~ -r± (4) 

a 3 

N s ~ r— (N s f s ) 2 (5) 



D 

.2 > T 1 o? 



k s~ m D-^i with —~—NJ s (6) 



>s IU D 

1 a 2 

f-NJ 

r fc 



Z - T^J, (7) 
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In (5)- (7) we have included the usual final state occupation factors in the 
approximation 1 + f s ~ f s . These occupation factors are essentially identical 
to those appearing in the early stages of the bottom- up picture [1] 
Eq.(4) says that the soft particles dominate the Debye mass calculation and 
that they do so in a self-consistent manner. Eq.(5) says that the soft parti- 
cles are produced from other soft particles (produced a little earlier) by the 
Bethe-Heitler process. Eq. (6) says that once particles are produced at a mo- 
mentum mo, by the Bethe-Heitler process, they then multiply scatter with 
soft particles produced a little earlier to reach the momentum k s given in (1). 
Finally, Eq.(7) says that a soft particle scatters with momentum transfer k s , 
with other soft particles once in the time period r. Thus, these particles are 
at the edge of equilibration. 

Now, since we have not been able to definitively determine that the mo 
in (1) is a screening mass rather than instability mass why have we assumed 
it to be a screening mass? The answer is straightforward; it is not internally 
consistent to suppose the m D of Eq.(l) to be an instability mass. To demon- 
strate this we need to calculate r m , the time over which particles produced 
at a scale vtid stay at that scale before having their momentum increased (by 
say a factor of 2), and f m the occupation number of particles at the scale m. 
Clearly 

1 a 2 Ar . 1 , , 

- N *f* = xi/fi+A/m ( 8 ) 



and 



The fact that r m ^> 1/tod while af m <C 1 means that the particles at mo 
do not come from an instability, because the instability would quickly raise 
af m to a value at least on the order of 1. 

Finally, what is the difficulty in deciding whether the particles described 
by Eq.(l) have a sufficiently symmetrical distribution so that screening dom- 
inates instability? The heart of the problem can be seen using Eq.(2) to 
calculate the "Debye mass" at r due to particles produced at r . One finds 

■2, , "N.{r,r a ) QUQ,T ) 3S " , 1n . 

m ° (T - To) ~n^r~ w .r) • (10) 
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Of course if t /t <C 1, m 2 D {r, r ) will be an instability mass, because the par- 
ticles produced at r will have a very asymmetric momentum distribution at 
r [3]. But then if tq/t <C l,m 2 D (r,To) <C iti 2 d (t) so that these gluons are not 
important. The gluons produced very close to the time r have equilibrating 
interactions among themselves and their distribution may well be quite sym- 
metric. What is difficult to decide is the role of gluons with, for example, 
tq = 1/2t. These gluons have begun to have an asymmetric distribution and 
their contribution to m 2 D is not much less, depending on 5, than the glu- 
ons produced at r. Thus, there are likely various contributions to m 2 D some 
screening and some favoring instability. Larger values of 5 favor screening 
because the gluons produced later are more emphasized, however, we have 
not been able to find a technical argument which decides whether screening 
or instability dominates for a particular 5. 



3 Matching onto the bottom-up picture 



As observed in Ref.[3] the early time (1 < Q s t < a~ 3 / 2 ) part of the bottom- 
up picture is not correct because what was assumed to be a Debye mass was 
in fact a scale characterizing an instability. The late time (a~ 5//2 < Q s t) part 
of the bottom-up picture should be consistent because the Debye mass was 
there determined from thermalized soft gluons. It is not immediately obvious 
whether at intermediate times (a~ 3 / 2 < Q s t < a~ 5 / 2 ) the bottom-up scenario 
is self-consistent. Here the interaction of the soft particles among themselves 
is of the same strength as the soft particles with the hard particles. In each 
case the time scale for an elastic scattering with momentum transfer on the 
order of k b s ottom " up ~ Q s \/~ct is the time r itself. We do not know how to 
answer the question of the stability of the intermediate time region of the 
bottom-up scenario although this issue is crucial for determining exactly how 
our present solution matches onto the bottom-up picture. 

Our present solution is in fact a family of solutions, depending on the 
parameter S. We now observe that the solution given in Eq.(l) agrees with 
the corresponding intermediate time quantities of the bottom-up scenario at 
a single time given by 



when < 5 < 1/3. For reference we recall that when a 3 / 2 < Q s r < a 5//2 




(11) 
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the bottom-up quantities are 



N s ~ \i/9 » fc s ~ VaQs, a/a niZFi m/9 > m £> 



(g s r)V2' a V4(g sT )V2' (g s r)V4- 

(12) 

Thus when r = r there should be a transition from our present solution, 
Eq.(l), to the bottom-up solution, Eq.(12). For r > f the bottom-up solution 
should be the correct one. We note that in case 5 = 1/3 the transition is 
to the beginning of the final stage of bottom-up, starting at Q s r ~ a~ 5 / 2 . 
Of course this transition to bottom-up can occur, for 5 < 1/3, only if the 
intermediate time part of bottom-up is internally consistent. If this is not 
the case then S must be greater or equal to 1/3. 

If 5 > 1/3 the solution (1) changes character at a time T\ given by 

Q s n ~ (l/a)5W (13) 

at which time f s ~ 1. In this case our solution goes into an evolution much 
like the final phase of bottom-up where the soft gluons are thermalized and 
the harder gluons feed energy into the soft thermalized system causing the 
temperature to rise with time until the whole system is thermalized. The 
difference now is that N s (t,t ) is larger than the number of hard particles 
so that as time grows larger and larger the gluons from N s (t, t ) at smaller 
and smaller tq disappear into the thermal bath, through branching as in 
bottom-up. The equations governing this evolution are 

de N s (t,t ) 

— ~ T A — ~ k s (r ) (14) 

GST GST T 

where 



k s (r ) ~ a 4 T 3 r 2 (15) 

Eq.(15) is the same as in bottom-up where now k s (r ) plays the role of 
kbr while Eq.(14) is the same as (14) and (15) of the first paper in Ref. [1] 
when one replaces N h of that reference by N s (t, t ). Using (15) to give the 
dependence of r on T(r) and r one finds, using (2), 

35-78,5 15-36,5 

T ~ Q s a s^-m (Q s r) 39*- 1( > . (16) 
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It is straightforward to verify that when 8 = 1/3, T ~ a 3 Q 2 s T as in Ref.[l] 
and also that T 4 ~ a (Q t) ' ^ ne en dpoint of the non-equilibrium evolution, 
when Q s t = a~ 13 / 5 with T = Q s a 2 ^ independent of the values of 5. Thus 
when 5 > 1/3 the matching onto bottom-up occurs only at the final time, 

How big could 5 possibly be? Certainly 5 < 10/21 because when 5 = 
j[,N s k s is equal to the energy density of the hard particles and this is an 
absolute limit as to the amount of energy carried by the soft gluons. 



4 What happens to the instability? 

So far we have said little about the instability except to comment that it 
should turn into something else in a very short time. Here we try to be a 
little more quantitative starting from the Abelianization picture of Arnold 
and Lenaghan[8] who suggest, and give evidence, that the instability should 
grow in one color direction. Let f\ be the gluon occupation in that Abelian 
direction at a momentum scale to, the instability mass. Of course when 
there are many gluons in a single color direction they may scatter elastically 
into other color directions through a lowest order collision term. Call the 
occupation number in these other color directions, but still at a momentum 
scale to, f 2 . Then 



f ~(™ 3 /i)^[/i(l + / 2 )]. (17) 

The first term on the right-hand side of (17) is just the number of gluons 
in the Abelian direction, the second factor is the cross section for scattering 
into the other color directions, and the third factor gives the remaining boson 
factors. Write (17) as 

' d ^~m{ah) 2 . (18) 



l + / 2 dr 

Now we suppose that /i(r) ~ e mr /i(0). This gives 

/ 2 (r) - [a/i(r)] 2 (19) 

so long as aj\ is small. The lowest order collision term which leads to (19) is 
only systematic so long as afi < 1. Thus, although (19) allows one to follow 
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the rise of f 2 to become a quantity of order one it is not possible to reliably 
follow f 2 into the region where it becomes much greater than 1, although the 
elastic collision terms themselves clearly cause f'2 to become as large as j\. 
Our proposed solution supposes that the instability is damped by a growth 
of soft gluons in all color directions and that there is then a transition to 
a solution of the form which we have presented. In this regard the recent 
numerical calculations by Dumitru and Nara[13] are interesting in that they 
show that the Abelian instability seems to run its course before there is a 
significant change in the hard particle distribution. It may well be that what 
is being seen is a damping of the instability due to soft gluon production of 
the type indicated in (17). 
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